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We construct the 1/S spin- wave expansion for double exchange ferromagnets at T=0. It is assumed that 
the value of Hund's rule coupling, Jh, is sufficiently large, resulting in a fully saturated, ferromagnetic 
half-metallic ground state. We evaluate corrections to the magnon dispersion law, and we also find that, in 
contrast to earlier statements in the literature, magnon-electron scattering does give rise to spin wave damp- 
ing. We analyse the momentum dependence of these quantities and discuss the experimental implications 
for colossal magnetoresistance compounds. 
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The phenomenon of colossal magnetoresistance 
(CMR), with its potential technological applications, has 
motivated an extensive experimental and theoretical re- 
search effort directed at the understanding of properties 
of doped manganese oxides In particular, consider- 
able attention has been paid to magnetic properties of 
these compounds. An ubiquitous feature shared by the 
numerous theoretical models of the CMR compounds is 
the presence of a strong ferromagnetic Hund's rule ex- 
change coupling, Jjj, between the spins of itinerant Cg 
electrons and the core spins of Mn ions. The kinetic 
energy of itinerant electrons is then minimized when the 
ionic spins are parallel to each other; this gives rise to the 
conduction-electron mediated double exchange ferromag- 
netism Q . The physics of double exchange is thus com- 
pletely different from that of both Heisenberg exchange 
and RKKY interaction (the latter corresponding to the 
case of small Jh)- 

Since the value of core spin is relatively large, S ~ 3/2, 
studying the effects of double exchange interaction in the 
large-S* limit should provide at least a qualitative descrip- 
tion of the low-temperature properties of the manganites. 
At the same time, the quantum nature of the core spins 
does affect the behaviour of the system in a profound 
way, as indicated by the variational results. At finite 5, 
the continuum of Stoner-like single-particle excitations 
with finite energies is present even at Jh — > oo [^j; also, 
the spin wave spectra obtained variationally show appre- 
ciable quantum-spin corrections . These conclusions 
are fully supported by the exact diagonalization studies 
of one-dimensional double exchange systems 1^-0 . 

It is therefore somewhat surprising that, apart from the 
leading-order calculations [p|-pT| (which do not account 
for any quantum corrections), the perturbative spin wave 



theory of double exchange magnets remains underdevel- 
oped. The objective of the present paper is to partially fill 
this gap. Our approach (unlike that of Ref. [Q) remains 
valid in the experimentally relevant case of large carrier 
densities. We will show that the subleading terms in 
the 1/S expansion, which originate from magnon-electron 
scattering, provide corrections to the magnon dispersion 
law, and also give rise to magnon damping. We will also 
see how the momentum dependence of these quantities is 
affected by the Fermi surface geometry. The experimen- 
tal data will be discussed briefly. 

We start with the usual double exchange Hamiltonian, 



Jh_ 
2S 



PA 



(1) 



* E-mail: golosov@nordita.dk 



Here Cja are the electron annihilation operators, Si are 
the operators of the core (localized) spins located at the 
sites of a square (or simple cubic) lattice, and the vector 
is composed of Pauli matrices. We assume, for sim- 
plicity, that only one conduction-electron orbital (with 
two possible values of spin projection, a [) is avail- 
able per each site i. Throughout the paper we use units 
in which hopping t, Ti, and the lattice spacing are all 
equal to unity, and we consider the T = case. 

There is presently no reason to doubt that for 5 > 3/2 
the ground state of the Hamiltonian (|l|) is ferromag- 
netic, at least for an infinite system in two or three di- 
mensions and with a finite number of electrons per site, 
X < 1 (this is corroborated by the variational calcula- 
tions 1^,^,^; see also Ref. Q for the ID case, and Ref. 
Q for finite-size systems). To leading order in 1/S, the 
electron spectrum in the ferromagnetic state is given by 
el"^ = e- ip Jh/'^, where for the tight-binding model of 

Eq. (|l|), = — Sa=i cos/ca, and d is the dimensionality 
of the lattice ||l^ . We will only consider the half-metallic 
case [|6l when the chemical potential lies below the bot- 
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torn of the upper band, ep < \ Jh — d + 0(1/5'), so that 
only spin-up electrons are present in the ground state. 

The magnon operators Oi are introduced by means of 
the Holstein - Primakoff transformation (including the 
subleading terms) of the operators 5*^, and the canonical 
transformation pTI Ti —>■ Ti' = exp{—U)Ti exp{U) with 
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(where N is the number of lattice sites and the sum- 
mation is over the first Brillouin zone) is then applied to 
the Hamiltonian, thereby accomplishing the transition to 
the "true magnons" . The resultant Hamiltonian can be 
written as an expansion in powers of S^^^^: 
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Here, the term Ti'i oc S ^1"^ has been eliminated by the 
canonical transformation, and 
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Here, the subscripts 1,2,.. stand for pi,p2,..., and 
means that the quasimomentum conservation law is en- 
forced. The spin wave energy cup is equal to the on-shell 
value of the real part of the magnon self energy (cf. Ref. 
pof ) . The leading-order term in self energy, which orig- 
inates from 7^2, is real and coincides with the earlier re- 
sults |§,nG,ni|, which in the Jh oo hmit reduce to 
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where n^- is the Fermi distribution function for the spin- 
up electrons. Using the explicit form of e^, one obtains a 
Heisenberg-like expression, w?'' = \E\{d+ep)/2dS . Here, 
E = J e^^d'^q/ (27r)'' is the total energy of electrons cal- 
culated with respect to the center of the lower band. 

Evaluation of the higher-order terms in Eq. (^ in- 
volves repeated commutations of the operator U (see Eq. 
(^) with the original Hamiltonian, Eq. (^. In order to 
calculate the leading correction to the energy of a sin- 
gle magnon, we will need to collect only the terms which 
contain neither c^^ nor ct^ and are quadratic in magnon 

operators. Such terms do not occur in Ti.'^, and we find 
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Naturally, both H2 and 7^4 have well-defined Jh 00 
limits. In what follows we will for simplicity restrict our- 
selves to the case of infinite Jh ; straightforward general- 
ization to the finite- J// case can always be accomplished 
with the help of Eqs. (|) and (|). 

The first quantum correction to the magnon self energy 
is proportional to S~'^ and includes the first-order per- 
turbative contribution from Ti'^ as well as second-order 
contribution T,2{0,P) from (corresponding to the first 
diagram in Fig. || below) JlTI I: 
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For the X ^ 1 case, this agrees with Ref. fl^ . 

For a 2D system with band filling values a; = 0.3 and 
X — 0.4, the quantity ujp^^ S"^ is plotted in Fig. |^ a (solid 
and dotted lines, respectively). The dashed lines rep- 
resent the corresponding Heisenberg-like fits, uj^'^ S'^ = 
2D'^^\2 + ep)S^, where [R 
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is the first quantum correction to spin stiffness. The de- 
viation of uj^'^ from uj^^"^ is not large, but it increases dra- 
matically for smaller x, in agreement with the numerical 
Q and variational results. 

The doping dependence of the spin stiffness, D = 
+ for 5 = 3/2 is shown in Fig. b (solid 

line) . The dashed line represents the leading-order term, 
= \E\/8S. We see that for S = 3/2 the value of D^^^ 
is not small numerically [l9| , implying that quantum cor- 
rections to the spin wave spectrum cannot be omitted in 
any quantitative treatment [|,|. Note that D{x) is not 
symmetric with respect to quarter filling, x = 0.5, re- 
flecting the loss of particle-hole symmetry at finite S. 

Spin wave damping is given by the on-shell value of 
the imaginary part of magnon self energy, and the first- 
order perturbation theory terms obviously do not con- 
tribute to it. It is also easy to sec that to second order in 
1/5, lmS2{'^piP) ~ 1111^2(0, p) = 0. This is because the 
energies of all the intermediate states occurring in the 
second-order perturbation theory terms are higher than 
that of a single-magnon state (which to leading order 
equals the ground state energy). Therefore, the integra- 
tion contour in the energy space terminates at the pole 
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(cf. Eq. (||) with w = 0), so that the latter does not give 
rise to any imaginary contribution. In the third order 
in l/S*, a multitude of new diagrams appear, combining 
all possible vertices from 7^2, H3, and H'^. By the above 
reasoning, their leading-order imaginary parts vanish on 
shell [^, with one exception. This exception is the sec- 
ond diagram in Fig. |2| which corresponds to the term 
with the second-order pole in the standard third-order 
perturbation theory formula |Q. This is nothing but a 
self energy correction to S2(tt'p-,p) (Eq. (||)), and the net 
third-order contribution to the magnon relaxation rate is 
thus given by the diagram on the r. h. s. in Fig. ^ [ p2| : 

x^(-f (10) 

Note that J(p, r) is of the order of 1/ S. 

Thus in a double exchange system, the spin wave 
linewidth remains finite even in a Jh — > 00, T = case. 
While such a possibility was discussed earlier in con- 
nection with recent experiments |p3| , it is at variance 
with the prevalent view |^,^. Magnon damping in a 
double exchange ferromagnet is due to magnon-electron 
scattering: a magnon can scatter into a lower-momentum 
state, while exciting an electron from under the Fermi 
level (see Fig. The momentum transfer into the elec- 
tronic subsystem allows us to speculate that an electrical 
current may arise in this process. 

The value V{p) for different values of bandfilling x in 
the 2D case is plotted in Fig. |l] c. It shows strong mo- 
mentum dependence, reaching (for S = 3/2) up to 10 % 
of the corresponding magnon energy. We note a differ- 
ence from both a Heisenberg ferromagnet (with r(p) = 
at T = 0) and an RKKY (s-/ exchange) ferromagnet, 
where the leading term in V{p) is threshold-lik e p5| . 

We note that Eqs. (p|-p^) (as well as Eqs. (g-|8|)) re- 
main unchanged if the first term in the Hamiltonian (|l|) is 
written as Vlr erd c,- with an arbitrary er. It is thus 

^k,a k ka„ka k 

possible to use Eqs. (|g|-[lC|) to evaluate r{p) for an arbi- 
trary single-band electron dispersion. The overall profile 
of r{p) (as well as that of cujj) is sensitive to both the 
carrier concentration value x and the details of electron 
bandstructure throughout the Brillouin zone. There are, 
however, two generic features in the behaviour of T{p): 
1. Long- wavelength limit. When p is small in com- 
parison to the Fermi momentum kp, T{p) is propor- 
tional to p^ in the 2D case and to p^ in three dimen- 
sions Here, the factor p^ originates from the ma- 
trix element {e^ — ip+qY Eq. (p^, and is multiplied 
by a p'^ corresponding to the space volume available to 
the virtual magnon with momentum p — f and energy 

remaining factor 



of p^ comes from electron kinematics. In the one dimen- 
tional case, T{p) = everywhere at p < kp. 
2. The anomaly at p — kp. The Fermi surface geome- 
try manifests itself through a weak singularity of r(p) at 
p — kp |2^. Indeed, the quantity J{p,r) (see Eq. (p^)) 
has a singularity at r = Ikp, and vanishes at io'^^lp — ■ 
The anomaly in T{p) is due to the tangency between these 
two surfaces (in r-space), which occurs aX p — kp. At this 
point, the second (in 2D) or third (in 3D) derivative of 
r(p) suffers either a jump or a logarithmic divergence, 
depending on the local geometry of these two tangent 
surfaces. These singularities are too weak to be visible in 
Fig. 1^ c, except for x — .49, when owing to the flatness 
of the Fermi surface the singularity acquires a nearly-lD 
(jump in the first derivative) character. 

The anomaly in T(p) is of course accompanied by a 
singularity in the third-order (in 1/S) term in up. It is 
contained in the real part of the same diagram, and for an 
isotropic dispersion law amounts to a jump in the second 
derivative (2D) or to a logarithmic divergence in either 
the first (ID) or third (3D) derivative of ojf. 

In view of the anticipated Fermi surface geometry of 
the CMR compounds |^^, we mention the case when 
the Fermi surface is nearly-cubic (nearly square). As 
the magnon isoenergetic surfaces are unlikely to have the 
same shape ||2^, this case is not reducible to ID. In fact, 
when p is perpendicular to a flat part of the Fermi sur- 
face, the singular term in r{p) is proportional to {p — 
kpf^eip-kp) in 2D and to {p-kpfe^p-kp) in 3D [|7)- 
The singularity in ojjj is given by —{p — kp)'^/^9{p — kp) 
in 2D and by {p — kp)'^\n\p — kp\ in 3D. 

Some comments are in order concerning the experi- 
mental situation. Low temperature spin dynamics, in- 
cluding spin wave damping and the deviation of magnon 
spectrum from the Heisenberg-like form, is studied in- 
tensively in both 3D perovskite and quasi-2D 
layered |^ CMR manganites. In both cases, unexpect- 
edly large low-T values of T{p) at large p were found, 
which is consistent with our results (in some compounds 
r(p) may be due in part to electron-phonon cou- 
pling). However, the lack of detailed knowledge of elec- 
tron bandstructure restricts one's options in comparing 
theoretical and experimental results quantitatively. 

The low- momentum measurements of T(p) in the per- 
ovskite manganites |2|] yield the p'^[TlTi{T/uip)]'^ depen- 
dence, attributable to the magnon-magnon scattering 
psf . We suggest that these measurements, both in 2D 
and 3D compounds, should be performed at lower tem- 
peratures in order to reveal the p^ or p^ contribution of 
magnon-electron scattering (which is likely to have a nu- 
merically small prefactor |^^). One should also try to 
identify the p ~ kp anomalies in F and uj (the irregu- 
larity seen in the data of Ref. at p ~ 7r/2 may be a 
possible candidate). We also suggest that a systematic 
study of doping dependence of spin stiffness should be 
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performed, in particular, in the layered compounds |35[ . 

In conclusion, we have calculated the spin wave 
linewidth and correction to the spin wave energy for 
a double exchange half-metallic ferromagnet. The pro- 
posed new measurements would verify to what extent the 
double exchange model accounts for the low-temperature 
spin dynamics of the CMR compounds. 
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Note added in proof: Our Eq. (0) agrees with Ref. @ , 
reporting a 1/S expansion for the model with a strong fer- 
romagnetic exchange, J/, between the core spins. How- 
ever, the second-order calculation in itself is not suf- 
ficient to evaluate the magnon damping at J/ = 0, given 
by Eq. (0) above. 
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FIG. 1. (a) The leading 1/5- correct ions to the spin wave 
energy in a 2D system with carrier concentration x = 0.3 
(solid line) and x = 0.4 (dotted line) as functions of momen- 
tum. The dashed lines represent respective Heisenberg-like 
fits, ili^^^S^. (h) Doping dependence of spin stiffness D in a 
2D system with S = 3/2. The dashed line corresponds to the 
classical value, D'^^^ . (c) Momentum dependence of spin wave 
damping, r(p), for a 2D system with x — 0.49, x — 0.4, and 
X — 0.3 (solid, dotted and dashed lines, respectively). For 
X = 0.49, note the p = kp anomaly, visible at p~ {''r/2, 7r/2} 
(see the inset), and the smallness of r(2fcF). 



FIG. 2. Magnon self energy diagrams. Dashed and 
solid lines are the unperturbed Green's functions of magnons, 
[uj + iO)~^ , and of the spin-up electrons, respectively. The 
bold dashed line is the exact magnon Green's function, which 
to required accuracy is given hy {lu — lu +iO)~'^. The boxes 
represent the interaction, i.e. the c|c|a^a-terms in Eq. (^. 
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